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Data on low-energy hadron production in collisions of nucleons, pions and kaons with light nuclei are needed 
for many astrophysical and accelerator applications. Modern simulations have reached a level of accuracy that the 
lack of detailed understanding of hadron production processes has become one of the most important limitations 
to further improvements. After giving some examples of hadroproduction processes in astrophysics and neutrino 
experiments we briefly review existing fixed-target data on light nuclei. Preliminary results and prospects of 
current CERN measurements (HARP, NA49) are discussed. 



Introduction 

The increase in computational power in recent 
decades has allowed some very precise experimen- 
tal measurements to be made. In particle physics, 
detector effects have been simulated using the 
Monte Carlo technique to the sub-percent level 
to obtain accurate corrections to measurements. 
Effects such as particle decay dynamics, electro- 
magnetic shower propagation and energy loss are 
understood to a sufficient accuracy to allow these 
precise corrections to be determined. The interac- 
tions of hadrons, on the other hand reveal a con- 
siderable gap between the accuracy available in 
simulation programs and what is desired for cer- 
tain applications. The reason for this is because 
the underlying theory of the strong interaction, 
Quantum Chromodynamics (QCD), is computa- 
tionally difficult except for processes with large 
momentum transfer where perturbation theory 
can be applied. 

The authors of simulation programs for hadron 
interactions have adopted either an empirical ap- 
proach (insert tables and parametrizations of ac- 
celerator data into the computer code, for exam- 
ple, |H2p or a semi-empirical approach (design 
models which are motivated from a theoretical 



*Based on a talk given at the XIII International Sym- 
posium on Very High Energy Cosmic Ray Interactions, 
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point of view and then use accelerator data to 
tune the free parameters in the model to make 
the simulation fit the data, for example, |3I4| 1. 
These procedures both suffer from a lack of avail- 
able hadron production data. As described be- 
low, good data is only available in a small frac- 
tion of the kinematically allowed and physically 
important region and for a small sample of target 
nuclei. Model builders are frequently faced with 
a choice in how to perform the extrapolation into 
uncharted parts of the phase space. 

This article begins by describing some appli- 
cations drawn from both cosmic ray physics and 
neutrino physics where more accurate simulations 
of hadron production would considerably improve 
the interpretation of existing data in terms of fun- 
damental parameters or observables. We continue 
with a review of the current state of hadron pro- 
duction measurements; the authors' brief was to 
cover CERN experiments, although other experi- 
ments are also mentioned. 

Muon production in extensive air showers 

Extensive air showers (EAS) are produced by 
high energy cosmic rays colliding with air nuclei 
in the upper atmosphere. The accurate simula- 
tion of the resulting cascades may be used to pre- 
dict the response (e.g. the distribution of muons 
on the Earth's surface). These simulations can 
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Figure 1. Energies and initiating particle type 
of interactions which lead to muon production in 
showers of E = 10 14 eV (top) and E = 10 16 eV 
(bottom). The upper horizontal axes show the 
equivalent hadron-nucleon cm. energy. 



then be used as a detector calibration; to deduce 
the characteristics of the original cosmic rays (e.g. 
the primary energy) based on the measurements 
in a detector. While the original cosmic ray might 
be of an energy which exceeds by far the energy 
available in accelerator based experiments, the 
showers contain a considerable number of lower 
energy interactions. The number of muons in 
EAS at detector level, often used together with 
the number of electrons to determine the primary 
energy and particle type of the shower, is partic- 



ularly sensitive to hadron production processes in 
the energy range from a few tens of GeV to sev- 
eral TeV 0. 

Figure n shows the energy distribution of those 
hadronic hadronic interactions in EAS that pro- 
duce observable muons. For each muon arriving 
at detector level (1032 g/cm 2 ) with an energy 
> 250 MeV, one entry is made at the energy 
of hadron /iqm 
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li + x', (i) 



where /igm and /im can be considered as the 
muon's grandmother and mother particles, re- 
spectively. The simulations are done with COR- 
SIKA using QGSJET as high-energy 
hadronic interaction model and GHEISHA for 
interactions below 80 GeV. The curves illustrate 
the relative importance of different hadron inter- 
actions at various energies. Due to the competi- 
tion between interaction and decay processes, the 
most probable energy is always between 100 to 
200 GeV for charged pions and somewhat higher 
for charged kaons, almost independent of the pri- 
mary EAS energy. This energy range of inter- 
actions relevant to low-energy muon production 
in EAS is in the reach of fixed-target accelera- 
tor experiments. Uncertainties in the prediction 
of muon numbers due to modeling of low-energy 
interactions are presented in Refs. |8I9| . The un- 
certainties increase with lateral distance from the 
shower core and exceed, for example, 20% at 2 km 
distance and 10 19 eV. Therefore, a considerable 
gain in the quality of EAS simulations can be 
expected with new hadron production measure- 
ments. 

Neutrino production in the atmosphere 

Atmospheric neutrinos have become an excit- 
ing study as the Earth, it turns out, is just about 
the most convenient size to see neutrino oscilla- 
tions |10| by comparing the zenith angle and en- 
ergy distributions with calculations (for a review, 
see ^J)- Even with an unsophisticated calcu- 
lation, it is clear that a deficit of upward going 
muon neutrinos exists which is interpreted as the 
oscillation of — > v T . 

Cosmic rays with energies extending down from 
EAS energies to a few GeV collide with the at- 
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Figure 2. Response for sub-GeV muon-like 
events in SuperKamiokande to the energy of the 
primary cosmic ray nucleons (from 52)- Shown 
are simulations for solar minimum and no geo- 
magnetic cutoff, events from lower hemisphere 
(upward going leptons), and events from upper 
hemisphere (downward going leptons). 



mosphere. Neutrinos are produced in the decay 
of muons, pions and kaons in the resulting cas- 
cades. The neutrinos which are observed as con- 
tained events in underground detectors [E v ~ 
1 GeV) are produced mostly by the decay of pi- 
ons from primaries in the range 2-200 GeV, see 
Fig [3 Higher energy neutrinos can be observed 
as upward going muons crossing the detector and 
are produced by primaries with energies up to 
~ 10 TeV. At these higher energies, kaon produc- 
tion becomes extremely important due to their 
shorter lifetime than pions (which makes them 
more likely to decay than interact). Additionally, 
the energy available in the center of mass is much 
larger for K — > \w than for 7r — > [iv and the 
steeply falling cosmic ray spectrum increases the 
importance of these kaon decays. With the recent 
improvements in primary flux measurements, the 
lack of comprehensive hadron production mea- 
surements over a more extensive fraction of the 
phase space is currently dominating the uncer- 



tainty of the atmospheric neutrino fluxes. 

Accelerator neutrino beams 

Accelerator neutrino experiments have been 
used to extensively study the interactions of neu- 
trinos including the discovery of neutral cur- 
rents, precision electroweak measurements and 
the search and study of oscillations. The acceler- 
ator beams are produced in a similar way to the 
atmospheric neutrinos described above. Protons 
from the accelerator hit a target producing pions. 
The neutrinos produced in the decay of the pions 
(mostly Vfj) are used for the experiment and the 
beam is stopped before many of the muons can 
decay. The neutrinos which are produced have 
a broad energy spectrum and are mostly how- 
ever there is a contamination of z/ e , T7^ and from 
muon and kaon decay. Precision measurements 
rely on an accurate knowledge of the neutrino 
energy spectrum and this is obtained from simu- 
lation with similar limitations from hadron pro- 
duction as for the applications described above. 
Many of the hadron production experiments have 
been performed by neutrino physicists with a spe- 
cific goal of understanding their neutrino beam. 
For a review, see Ref. |13| . 

A radical new way of producing very intense ac- 
celerator neutrinos has been proposed (a neutrino 
factory) in which a beam is first produced in the 
way described above from a low (~ 2-10 GeV) 
energy proton beam. After the pions have been 
allowed to decay, the muons are collected and 
then inserted into a further accelerator where 
they are accelerated to energies ~ 50 GeV. 
These muons are then circulated in a storage ring 
with long straight sections pointing towards sev- 
eral detector sites. The decay of the muons in 
the straight sections produces an intense source 
of neutrinos of known composition and spectrum. 
The gain in the concept is that power is applied 
by the accelerator directly to the muons. There 
are many difficulties to overcome, such as the con- 
struction of the target region, and the method of 
inserting the muons into the accelerator. It is pos- 
sible in the further future that the concept may 
permit a muon collider to be constructed. One 
of the purposes of the HARP experiment, to be 
described below, is to investigate the hadron pro- 
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duction in order to optimize the efficiency of a 
neutrino factory target. 

To summarize, there is considerable interest in 
the improvement of hadron production models for 
many applications. 

Currently available data 
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Figure 3. Summary of existing measurements of 
hadron production as a function of primary mo- 
mentum pprim and secondary pion momentum p v . 
The black (red) boxes indicate the importance of 
each phase space region for downward (upward) 
atmospheric neutrino flux at SuperKamiokande. 
The difference is due to the Earth's magnetic 
field. The vertical bars indicate where experi- 
mental data exist on low- A nuclei. The indi- 
cated measurements are Abbott et al. Cho 
et al. HS|, Eichten et al. [THj, Allaby et al. Q7|, 
Barton et al. JS], Serpukhov [T§|, Atherton et 
al. HJ, and SPY |3T]. 



Figure [3] shows the current status of hadron 
production measurements for proton projectiles 



in the region below 1 TeV primary energy. The 
vertical lines indicate the locations as a function 
of primary energy E p and secondary energy E 
where some measurements have been made. An 
additional degree of freedom, not indicated on the 
plot is the transverse momentum px of the sec- 
ondary particles (the production is symmetric in 
the ip angle). For most of the applications de- 
scribed above, the total production as a function 
of E p and E is of main interest (i.e. integrated 
over pt ) . Given the limited number of data points 
this integration is ambiguous and depends on the 
functional form used for the extrapolation to un- 
measured phase space regions (for example, see 
discussion in |22|). 

Since the pr distributions control the lateral 
extent of a shower, the functional form of the 
hadron yields as a function of px is also impor- 
tant for EAS simulations. Atmospheric neutrinos 
are only observed one at a time (i.e. the chance 
of detecting two neutrinos from the same cosmic 
ray are negligible), the form of the pt distribu- 
tion is therefore less important for atmospheric 
neutrinos, however it is important to have the 
yield integrated over the full region of px- Neu- 
trino beams select and focus pions from specific 
regions of phase space and so px distributions are 
important here also. 

Many of the measurements shown in figure |3 
have been made using 'single armed spectrom- 
eters'. These are devices where the secondary 
particles are measured with a string of magnets 
tuned to transmit particles of a given momentum 
at a time. Examples of experiments of this type 
are given in Refs. |lfil2QI2i| . The experiments 
can be 'ready made' by using the magnets in a 
secondary beam at an accelerator center such as 
CERN as the measuring device. A big advantage 
of this technique is that precise particle identifi- 
cation tuned to the momentum which the spec- 
trometer is being used can be made to separate n, 
K and p extremely well. The measurements are 
somewhat painstaking, since only one angle and 
momentum of the secondary particles can be mea- 
sured at a time and this is the reason why these 
measurements do not yet cover sufficient points 
in phase space for building models of hadron pro- 
duction. 
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A second technique is to use an emulsion stack 
in an accelerator beam. These experiments are 
limited in statistics and particle identification, 
however they can measure the complete phase 
space region in one exposure. Bubble chamber 
experiments have also been used, e.g. |23| . 

The experiments described below aim to mea- 
sure a large fraction of (or the complete) sec- 
ondary phase space in one setting by using large 
acceptance detectors. This will allow complete pt 
distributions to be obtained in large regions of E. 
The experiments arc not quite so good at particle 
identification as the single- arm experiments. 

HARP 




beam 



Figure 4. Layout of the HARP experiment |241 

Eg. 



The HARP hadron production experiment [21] 
was built specifically to address the lack of cover- 
age of phase space with good hadron production 
data for the applications as highlighted at the be- 
ginning of this article. A comprehensive range of 
targets (LH 2 , LD 2 , Be, C, LN 2 , L0 2 , Al, Cu, Sn, 
Ta, Pb, H 2 0) of different atomic weights A and 
thickness have been exposed to beams from the 
CERN proton synchrotron (PS) with momenta 
between 1.5 and 15 GeV. Beams of both posi- 
tive and negative targets have been used. Several 
replicas of neutrino beam targets have also been 
exposed. There are several solid targets with A 



close to that of air and cryogenic (liquid) nitrogen 
and oxygen targets have also been exposed. 

The beam is produced from a primary PS tar- 
get at 24 GeV and the secondary particles are mo- 
mentum selected and tagged with gas Cherenkov 
and time of flight detectors in the beam. All inter- 
actions are triggered and the type of beam parti- 
cle determined offline; thereby, interactions with 
parent protons, kaons and pions are all available. 
The beam instrumentation also measures the tra- 
jectory of the incoming particle so that the impact 
position on the target can be reconstructed. 

The detector is shown in figure 0] There are 
two distinct groups of detector to select different 
regions of phase space; large angle and small angle 
with respect to the beam direction. 

The large angle detectors consist of a 80 cm di- 
ameter by 1.5 m long time projection chamber 
(TPC) enclosed in a 0.5 T solenoidal magnetic 
field. A hole in the center of the TPC allows the 
target to be inserted in the center of the sensitive 
region of the detector so that backward going par- 
ticles can be measured. The TPC is the primary 
detection device for large angle tracks and pulse 
height information from the detector allows parti- 
cle identification through ionization loss (dE/dx) 
measurement J2B]. A layer of resistive plate detec- 
tors (RPCs) with ~ 160 ps resolution surrounds 
the detector for particle identification by time of 
flight The trigger is provided by a series 

of scintillating fiber detectors mounted along the 
central hole of the TPC. 

The tracking detectors in the forward region 
(small angle) are drift chambers with three views 
(vertical and ±5°). A vertical magnetic field 
bends the particles in the horizontal plane to de- 
termine the particle momentum. Downstream, a 
plane of scintillation counters with ~ 200 ps tim- 
ing resolution is used to identify the particles by 
time of flight and provides 3<r ir/p separation up 
to 4.5 GeV. A large (31m 3 ) Cherenkov detector 
filled with C4F10 provides good 7r/p separation at 
high momentum. Electron identification is pro- 
vided by a scintillating fiber and lead sampling 
calorimeter. 

The different phase space regions with the 
corresponding detector components employed for 
particle identification are shown in Fig. [5] The 
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Figure 5. Secondary particle distribution in p- 
Be interactions at 15 GeV. Phase space regions 
for different methods of particle identification are 
superimposed. 
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squares show the expected secondary particle dis- 
tribution for 15 GeV protons on Be. 

Preliminary data on a thin 20 mm (5% of an 
interaction length) aluminium target, a short sec- 
tion of a replica target for the K2K experiment 
have become available recently I^Hj. Figure 
shows the yield of secondary particles identified 
as pions as a function of both the momentum p 
and opening angle 9 Z in the horizontal (bending) 
plane for momenta above 0.2 GeV/c. Fiducial 
cuts were applied at 25 < \6 Z \ < 200 mrad and 
\9 y \ < 50 mrad, where 9 y is the opening angle in 
the vertical plane. 

NA49 

The NA49 experiment was first commis- 
sioned in 1994 and has been used since then with 
both lead ion and proton beams from the CERN 
SPS. The experiment is designed to cope with the 
extremely high multiplicity environment of a Pb- 
Pb collision at 400 GeV/Z. 

NA49 is installed in the H2 beam line at the 



Figure 6. Recent HARP hadron production data 
on a thin aluminium target (K2K target replica) 
as a function of momentum and of angle in the 
horizontal plane |28j. 



CERN SPS. The experiment (Fig. 0) consists 
of two (2mx2.5mx0.98m overall dimensions) 
TPCs (VTPC1 and VTPC2) supported inside 
1.7 T super-conducting magnets. Downstream, 
there are two further TPCs (3.9mx3.9mxl.8m) 
(MTPCL and MTPCR) which measure the tracks 
once they emerge from the magnetic field. Ap- 
proximately 100 samplings of the ionization loss 
are made per track in the combined TPCs to de- 
termine the particle identification from the rela- 
tivistic rise section of the energy loss curves. This 
is the principle form of particle identification in 
the experiment. The target is mounted in front 
of the first TPC, so, unlike HARP, only forward 
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Figure 7. Plan layout of the NA49 experiment. The four TPCs are shown (VTPC1,2 with magnetic field 
and MTPC L,R). The position of the target is indicated by a 'T' in the inset a) immediately in front of 
VTPC1. Beam position detectors (BPD) and trigger counters SI, 2, 3 are also shown. 



tracks may be measured. The detector has excel- 
lent acceptance for particles with Xf > 0. Parti- 
cles with momenta above 80 GeV/c, pass through 
the gap between TPCs (essential for heavy ion 
running to avoid space charge from ionization 
from the uninteracted beam). A new thin gap- 
TPC was added to the experiment to allow track 
detection (but not particle ID) in this region. 

A differential beam Cherenkov detector (two 
compartments, with pressures set appropriately 
to be just above and below the Cherenkov thresh- 
old of the desired particle) is used in the trigger 
so that only events which originated from protons 
are in the data sample. A beam particle tracking 
system (BPD) allows the impact position on the 
target to be known precisely. 

The experiment has recently been run with a 
low-bias trigger with a carbon target for the pur- 
pose of measuring hadron production for atmo- 
spheric neutrinos. A one week run with protons 
at 158 GeV and 100 GeV has yielded a sample of 
data with over 500,000 tracks [3D] and the analy- 
sis will be completed shortly. 

The NA49 experiment has finished its data tak- 
ing programme. Currently a new collaboration 
is forming to perform an extended series of mea- 



surements , including hadron production mea- 
surements at CERN. The experiment will use the 
existing NA49 apparatus and it is intended to col- 
lect a considerably more extensive set of measure- 
ments than obtained in the pilot experiment de- 
scribed above. It is planned to perform a primary 
energy scan from 40 to about 200 GeV and a vari- 
ety of projectiles will be studied including protons 
and light nuclei such as helium and carbon. 

Summary 

There are many reasons why the world's hadron 
production data is not enough. In particular se- 
ries of many points in phase space are required 
to allow the functional form of hadron produc- 
tion measurements to be known precisely across 
all pt and secondary particle momenta. Two 
new measurements at CERN, HARP (a purpose 
built experiment) and NA49 have been described. 
Data at both experiments have been collected and 
are currently being analyzed. Preliminary results 
from HARP have recently become available. 

Additionally, there are new experiments in 
the USA which have, or soon will have hadron 
production data. The E910 experiment at 
Brookhaven has data in the range 6-18 GeV 
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on beryllium |33I34| and the MIPP experiment 
(specifically built for hadron production measure- 
ments) at FNAL [35] has concluded its engineer- 
ing run in 2004 and is scheduled to collect data 
in the primary energy range 5-120 GeV in 2005 

M- 
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